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Abstract 

It is shown that for thin metalhc films thickness of which does not exceed thick- 
ness of skin layer, the problem allows analytical solution. In the field of resonant 
frequencies the analysis of dependence of coefficients of transmission, reflection and 
absorbtion on an electromagnetic wave is carried out. Dependence on pitch angle, 
thickness of the layer and coefficient of specular reflection and on effective elec- 
tron collision frequency is carried out. The formula for contactless determination 
(calculation) of a thickness of a film by observable resonant frequencies is deduced. 
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PACS numbers: 73.50.-h Electronic transport phenomena in thin films, 73.50.Mx 
High-frequency effects; plasma effects, 73.61.-r Electrical properties of specific thin 
films, 73.63.-b Electronic transport in nanoscale materials and structures 

1. Introduction 

The problem of interaction of an electromagnetic wave with the metal 
film attracts attention to itself for the long time already |I] - It is 
connected as with theoretical interest to this problem, and with numerous 
practical appendices as well [6] and jTj. 
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Researches of interaction of an electromagnetic wave with metaUic films 
were carried out basically for a case of specular dissipation of electrons 
on a film surface. It is connected with the fact that for more general 
boundary conditions the problem becomes essentially complicated and does 
not suppose the analytical solution generally. 

In the present work it is shown that for thin films, a thickness of which 
does not exceed a thickness of a skin - layer, the problem allows the an- 
alytical solution. In previous work [8] the case when the frequency of 
electromagnetic wave is less than plasma (Langmuir) frequency was con- 
sidered. 

Let us note, that the most part of reasonings carried out below would be 
true for the more general case of conducting (in particular, semi-conductor) 
film. 

2. Problem Statement 

We consider the thin slab of conducting medium on which the electro- 
magnetic wave falls. We denote the pitch angle by 9. We will assume, 
that the vector of magnetic field of electromagnetic wave is parallel to the 
surface of the slab. Such wave is called j9-wave (see [3] or [1]). 

We take Cartesian coordinate system with origin of coordinates on one 
of the surfaces of a slab, with axis X. directed deep into the slab. We direct 
the axis y parallel with the vector of magnetic field of electromagnetic wave. 
Under such choice of the system of coordinates the electric field vector and 
magnetic field vector have the following structure 

E = {E^{x, y, z), 0, E,Xx, y, z)} and H = {0, Hy{x, y, z), 0}. 

Further components of electric and magnetic field vectors we search in 
the form 

E,{x, z, t) = E,{x)e-'^'^'^''^'\ E^x, z, t) = E,{x)e-'^'+'^''^'\ 
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and 



-iujt+ik sin Oz 



Hy(x, z, t) = Hy{x)e 

Now behaviour of electric and magnetic fields of the wave in the slab is 
described by the following system of differential equations [3] 



dx 



ik sin 6 Ex + ikHy = 



47r 

< ikEx — ik sin 9Hy = — j^; 

c 



(1) 



dHy ^ 47r . 
— + ikE^ = — j^. 
dx c 

Here c is the velocity of hght, j is the current density, k is the wave 
number. 

We denote the thickness of the slab by d. 

The coefficients of transmission T, reflection R and absorption A of the 
electromagnetic wave by the slab are described by the following expressions 

ra, m 

1 .o 

(2a) 



7^^£|p(l)_p(2)|2 
4I I ' 



and 



41 I 



A=l-T-R. 



(26) 



(2c) 



The quantities P*--^^ {j = 1, 2) are defined by the following expressions 

cos^ + Z(^') 



J = 1,2. 



(3) 



cos^-ZW 

The quantity Z^^^ corresponds to impedance on the lower surface of slab 
under symmetrical configuration of the external magnetic field. This is the 
case 1 when 



Hy{0) = Hy{d), P,(0) = EM. EM = -EM 
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The quantity Z*^^-* corresponds to the impedance on the lower surface of 
slab under antysymmetrical configuration of the external magnetic field. 
This is the case 2 when 

Hy{0) = -Hy{d), E,{0) = -E,{d), E,{0) = E,{d). 

The impedance is thus defined as follows 

'"'-^y .-1.2. (4) 

We will consider the case when the width of the slab d is less than the 
depth of the skin - layer 5. Let's note, that depth of the skin - layer 
depends essentially on frequency of radiation, monotonously decreasing in 
process of growth of the last. The value 5 possesses the minimum value in 
so-called infrared case [ID] 

% = — , 

where Up is the plasma frequency. 
For typical metals [ID] 100 nm. 



Hence for the films thickness of which d is less than 6o our assumption 
is true for any frequencies. 

3. Problem solution 

The quantities Hy and Ez change a little on distances smaller than the 
depth of a skin - layer. 

Therefore under fulfilment of the given assumption d < 5 these electrical 
and magnetical fields will change a little in the slab. 

In case 1 when Hy{Q) = Hy{d) it is possible to accept that the value Hy 
is constant within the slab. 
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Variation of the quantity of ^/-projection of electric field on thickness of 
slab can be defined from the first equation of the system (1) 

d 

EM-EM-iMH,^ik^.ejE.a.. (5) 



From the second equation of the system (1) it follows that on the bound- 
ary of the film the following relationship is satisfied 

EM = EM = Hy^^^e. (6) 

The integral from the relation (5) is proportional to the value of the 
quantity of normal to the surface component of electrical field on the sur- 
face and therefore according to the relation (6) it's proportional to the 
quantity Hy. 

We define the coefficient of proportionality as 

d d 

G = — ^ / E,{x) dx = }. [ E,{x) dx. (7) 


With the help of (7) we rewrite the relation (5) in the following form 

E,{d) - E,{0) = {-ikd + ikGdsin^ e)Hy. 

Considering the antisymmetric character of the projection of electric 
field Ey in this case we receive 

Eyifd) = ikd{l - Gsin^e)^. (8) 

Therefore for the impedance we have 

zW = ^(l-Gsin2^). (9) 

For the case 2 when £^^(0) = Ez{d), it is possible to assume that z - 
projection of electric field Ez is constant in the slab. 
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Then the magnetic field change on the width of a slab can be determined 
from the third equation of the system (1) 

d 

47r f 

Hy{d) - Hy{0) = -ikdE, + — / jzix)dx. (10) 



Thus 

jz{x) = (t{x)Ez, 

where g{x) is the conductance that in general case depends on coordinate 

X. 

Let's introduce the longitudinal conductivity averaged by thickness of 
the slab, 

d d 

(Td = J jz{x)dx = ^J (y{x)dx. (11) 



Then the relation (10) can be rewritten with help (11) in the following 
form 

Hy{d) - Hy{0) = -ikdE, + ^^^E,. 
Considering symmetry of the magnetic field, from here we have 

HyiO) = \ikdE, - ^E,. 

For the impedance (4) we have 



= ?^ . (12) 

ickd — Airadd 

4. Coefficients of transmission, reflection and absorption 



From here according to (3) we receive expressions for the quantities 
P(^) (j = l,2) 
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2 cos^ + zA;(i(l — Gsin^ ^) 

p(2) ^ (47r(Jd - ikc)d cos e - 2c 
{4:7Tad- ike) d cos 9 + 2c' 
We will assume that length of a wave of incident radiation surpasses 

essentially the thickness of the slab, i.e. kd <^ 1. Then expressions (9) 

and (12) for impedances and expression (13) for quantity P(^) (j = 1, 2) 

become a little simpler 

Zi,^^ = --ikGd sin^ 9, = (14) 

° 2 '0 27Tadd ^ ' 

Substituting (14) into (3), we have: 

^(1) 2cos9 — ikGdsiv? 9 ^{2) 27Taddcos9 — c 
z cos 9 + ikGdsm'^9' ^ 27rM cos 6* + c 

Thus quantities R, T, A can be found according to the formulas (2). 

In a limiting case of non-conducting slab, when cr^f — > 0, G — > 1 from 

these expressions we have 

p(i) ^ _p(2) ^ '2 + ikd cos 9 
2 — ikd cos 9 ' 

from which 

r = i, R = 0, A = 0. 

Under almost tangent incidence when 9 — > 7r/2 we receive P*^^^ 
-1, P(2) ^ _i. Thus, we obtain that T ^ 0, P ^ 1, ^ ^ 0. 

Let the relation /c/ <C 1 be true. Then in a low-frequency case, when 
cj — > 0, the quantity cr^ for a metal film can be presented in the following 
form [TT] 



where 



1 3 X /" / 1 1 A 1 - e"""' , 



^{w) w 2w'^ J t^Jl—pe 
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Here / is the mean free path of electrons, p is the coefficient of specular 
reflection, ctq = aj^r/(47r) is the static conductivity of a volume pattern, 
T = l/vp is the time of mean free path of electrons, vp is the Fermi velocity. 

In a low-frequency case when the formula (16) is applicable, the co- 
efficients T, i?, A do not depend on frequency of the incoming radiation 
according to the formulas (2). 

For arbitrary frequencies these expressions are true under condition, 
that it IS necessary to use the rollowmg expression I — > - 



as a 



quantity I and the expression ctq — > 



ICJT 



instead of ctq. 



1 — liOT 

For the case kl <^ 1 the quantity G can be calculated from the problem of 
behaviour of a plasma slab in variable electric field, which is perpendicular 
to the surface of slab [il3j . 

Let's calculate the coefficients of transmission and reflection in the case 
when kd <^ 1. We will substitute expressions P(i) and P(2)^ defined ac- 
cording to (15), into the formulas (2). We receive 

^ Ai^"- r< a ■^'^ 2 



T = cos^ e 



^k^G sin' 6^-^ 
2 



c 



R = 



(cos e + ik^G sin' e){l + COS 9) 

ik-G sm^ 6 cos^ 6 

2 c 



(17) 



(cos 9 + ik'^G sin' 9){l + cos 9) 



(18) 



where the quantity G may be found from the solution of the problem of 
plasma oscillations |I21 and 

d 

(19) 



and eix) is the electric field. 



a 

G = — J e{x)dx, 



5. Coefficients of transmission, reflection and absorption in the 
case of specular reflection of electrons 
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In the case of mirror reflection of electrons according to (16) we have 
c"d = ctq, and the quantity G under the formula (19) can be calculated 



precisely, using the electric fleld in plasma layer constructed in |[T3 



At the proof of decomposition of the decision of an initial boundary 
problem in [13] the electric fleld in a metal layer has actually been con- 
structed 

q{x) - — + 2Aiy/ocosh(2:ox/y/o) ^ 
Aoo {ac - t]1)\'{t]q) cosh(zo/T/o) 
1 

2 J X+{r])X-{r])cosh{zo/r]) 

-1 

In this formula X{z) is the dispersion function from the problem of 
plasma oscillations, 

1 

A[z) = c + — j -T, ^ar, r/^ = ac, a = 



Z J — Z'^ VpX 
-1 

9aQ 2 S-Up aQiv — ico) 

ro is the Debaye radius, Ai = A(?7i) = — ac, ijI are the zeroes of the 
functions cosh.{zo / ij) , 

_ 2zoi _ 2ao{uj + zu) ^ - ±1 ± 2 • • • 

In the layer < x < d the electric fleld has the following form 

cosh[zo{2x - d)/r]id] Ai ''^^ r]f cosh[zo{2x - d) / ijld] 
cosh(zo/r/i) Vo X{r]l){r]f - r^) si^Hz^hl) ' 

The quantity G can be found easily by means of the equalities (19) and 
(20) (see also |I3]) and has the following form: 



Ai ^ 2Aiy/gtanh(zo/y/o) ^ Ai f tanh(zo/T/)T/^ dr] 
Aoo ZQ{ac- r]l)\'{T]Q) 2zq J A+(r/)A-(r/) 
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Here 



A (/i) = A(/i) ± ^-/i(r/i - /i' 



Integral from (21) we will calculate by means of methods of contour 
integration. Let us take advantage further of obvious equality 



1 



1 



1 



1 



Therefore the integral from (21) is equal to 
1 1 



A+(/i) A+(/i) 



1 f tanh(zo/r)r^(ir _ 1 

2 I 



-1 



A+(r)A-(r 



27rz 



1 



1 



A+(r) A-(r) 



tanh(2;o/r)r^ dr 



^2 _ ^2 



We take now a circle {'Jr : \z\ = R) with so big radius i?, that all 
finite singular points of function 

f( \ _ 'tanh(2;o/2;)2;2 

lay inside jr. Such points are points z = ±7/1, zeroes of the dispersion 
function X{z) are points z = ±770 (if (7, e) G ^^) (see [13]), and also 
polar singularity of the function tanh(2;o/2;). The last points are points 
z = r]l /c = 0,±l,±2,---. 

According to the theorem of the full sum of residues we receive 



fc=+oo 



Res + Res + Res + Res + Res + Res 



fc=- 



tanh(2;o/2;)^' 



27ri 



We note that 



X+(t) X-(t] 



^ tanh(2;o/^)^^ 
2=00 (^2 _ rjl)\{z] 



z2 - rH)\{z] 
tanh(2;o/T)r2 dr 



— r]\ 



Res 



tanh(zo/;^)^' 



:^^-r/?)A(^) ^oA(r/*)(r/f -r/?)' 
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tanh(£(j/-)-^ _ r/itanh(£()/77i^ 



tSLnh{zQ / z) z'^ ?7q tanh(zo/?7o) 
±.0 {z^ - r^l)\{z) = \'{mM-rf,) 



Res 



Thus, the integral from (21) is equal to 
1 

1 [ tanh(zo/r)r^ <ir Zq r]iiaii\i{zQ / rji) 



+ 



2 J A+(r)A-(r) Aoo 2\{r^i) 

-1 

277^ tanh(2;o/7yo) ^ 



2E 



Hence, the quantity G is equal to 



*4 



G = ^tanh^-:^ V (22) 



Let us note, that according to the evenness of the expression standing 
under a sign of the sum in (22), it is possible to simplify this sum and 
present the expression (22) in the following form 

r> \ k=+oo ^4 

Let us note, that the sum of the series from (23) is well approximated 
by the first member, i.e. instead of (23) it is possible to take 

Gi = — tanh ^ ^T- 

^0 m 2;2A(7yo*)(C-ryf) 

Let us carry out numerical calculations. We will enter the relative error 



Oi(n,£,(i) = 



G — Gi 



100%, Q = —,£ ^ 



UJp LUp 



G 

Then for a film of sodium (ojp = 6.5 • 10^^ sec~^, vp = 8.52 • 10^ cm/sec) 
of the thickness of 1 nanometer, 5 nanometers and 10 nanometers under 

oj = ojp and z/ = lO^'^cjp sec"^ accordingly we have: Oi = 1.42%, 1.38% and 
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1.98%. For G we have replaced an infinite series with the finite sum for 
= 10^ members. 

The quantity G is approximated by first two members of the decompo- 
sition (21) even more effectively, i.e. when we replace G with the quantity 

Ai ^ 2Air/g tanh(2:o/r/o ) 
Aoo z^{ac - t/o)A'(t/o 

The formula (24) means, that we have replaced the electric field (20) by 
two first components of Drude and Debaye, corresponding to the discrete 
spectrum. 

For the calculation of the quantity G^ explicit expression of zero of 
dispersion function 770 = ?7o(r^, e) is required. We will write the factorization 
formula of dispersion function (see [IB]) without the proof 

A(z) = Aoo(t/o' - z^)X{z)X{-z). (25) 

In (25) we introduce the following notations 

Aoo = A(oo) = ^ + ac - = ^(1 - - ieO), 

X(z) = -expy{z), V(z) = — y —-—dr, 



G{t) = Pf\, A±(r) = c'-ac- (r^ - ac)A±(r), 
A r 



Aj(r) = Ao(r) ± -rz, Ao(r) = 1 + - / — — = 1 + -In 



1 

TT , . . ^ T f dr' ^ r , 1 — r 



2' 2J t'-t 2 1 + r 

-1 



If we calculate values of the left and right parts of the equation (25) in 
the point z = then for square of zero of dispersion function after some 
transformations we receive following expression 



2 . A(z) ^ X{i) 



V{i)-V{-i) 
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TT 



Considering, that Ao(z) = 1 — — , we have 



A(z) = c2-ac+(l-^)(l + ac) = -^(^^ + 16^)+ (l-j^ [l + ^ie^-ie^) 



It is possible to present the function X(z) = -e^^^^ in the following form 



Let us find the sum 



1 f\nG{T)dT 



T — Z 



1 }\nG{T)dT ^ 1 }\YiG{T)dT 



— I 



2m 



1 }\nG{T)dT 1 }T\nG[T)dT 



27ri J T — i 

-1 



27ri 



f T\nG{T)i 
J r2 + l 



By means of these formulas we will transform the formula for the square 
of zero of the dispersion functions to the form 



or 



where 



2 , 2A(z) 
77o = -1 + — ; — exp 



A 



00 



1 }TGl{T)dT 

2^1 J r2 + l 



-1 



2 2A(z) 
n = -l + ^-^exp 

Ar 



*0O 



i rGi{T)dT 
nj r2 + l 





TT 



(3r2 - £2 + z£Q)(Ao(r) + -rz) + + 
G'/(r)= In 2 



TT 



(3r2 - £2 + ie^7)(Ao(r) - -ri) + ^^2 + 



Now a relative error 



G — 



100% 



(26) 
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for sodium films of the thickness of 1 nanometer, 5 nanometers and 10 
nanometers at cj = cu^ and v = lO^^cjp accordingly equals: O2 = 0.575%, 0.003% 
and 0.0004%. 

The graph of the relative error 02{^,£,d) as a function of the variable 
O for a film with the thickness of 10 nm at z/ = O.OOlcjpSec"^ is represented 
on the Fig. 1. 




Figure 1: Influence of continuous spectrum, d—lQ nm, v — 0.001a;p, 6 — 75°. 

The graph on Fig. 1 shows, that in the area oj cOp the contribution to 
the electric field, corresponding to the continuous spectrum, is insignificant, 
and can be neglected. Thus, function G is approximated by two composed 
Drude and Debaye, corresponding to discrete spectrum, according to (24). 

Using the formulas (17) and (18), with the use of the expression (24) 
for functions G, we will carry out graphic research of the coefficients of 
transmission, reflecion and absorption. 

6. Discussion of the results 

Let us consider the case of thin sodium film. We will construct graphics 
of the dependences of transmission, reflecion and absorption on quantity 
n = cj/cjp at pitch angle 9 = 75° (Figs. 2-10). 
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We will note, that near to a plasma resonance {cu ~ cUp) the coefficient 
of transmission! has a minimum, and reflecion and absorption coefficients 
have minimum values. Under a thickness of the film of 1.5 nanometers and 
at 1/ = O.Obujp in the areas of resonant frequencies {uj > cjp) all coefficients 
have one maximum more. Under increase of the thickness of the film from 
1.5 to 10 nanometers the second maximum vanishes. 

For a film with the thickness of 5 nanometers and under u = 0.02ujp in 
the field of resonant frequencies behaviour of all coefficients has so-called 
"edge" character ("paling"). Under the further increase of the thickness 
of the film frequency of "combs" increases also, and we can see growth of 
quantity of its "teeths" (Figs. 5 - 10). If for a film with the thickness of 
5 nanometers the quantity e = ly/ujp decreases, the quantity of "teeths" 
of combs grows sharply. 
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0.5 1.0 1.5 

co/co 

P 



Figure 2: Transmittancc, 6 = 75°. Curves 1, 2, 3 correspond to values of parameters d = 2 
nm, u = O.ObcUp; d = 5 nm, u = 0.030;^; d = 10 nm, u = O.ObcUp. 




Figure 3: Reflectance, 9 — 75°. Curves 1, 2, 3 correspond to values of parameters d — 2 
nm, u — 0.05a;p; d — 5 nm, u — 0.030;^; d — 10 nm, u — 0.05a;p. 




0.5 , 1.0 1.5 

CO/CO 

P 



Figure 4: Absorptance, 9 — 75°. Curves 1, 2, 3 correspond to values of parameters d — 2 
nm, 1/ — 0.05a;p sec~^; d — 5 nm, u — 0.03a'p, d — 10 nm, i/ — 0.05a'p. 




Figure 7: Absorptance, d — 5 nm, u — O.OOlcUp, 6 ~ 75°. 




Figure 10: Absorptance, = 10 nm, u = O.OOla;^, = 75°. 
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At the further increase of the thickness of a film frequency of teeths of 
the comb (number of linlcs of paling) increases. Let us note, that Fig. 10 
actually coincides with Fig. 3 of work [1], and Fig. 9 coincides with Fig. 
2 of [U. 

On Figs. 11 and 12 dependences of coefficients of transmission, reflection 
and absorption from the pitch angle 9 of the electromagnetic waves on a film 
are presented. These graphics show monotonously decreasing character of 
the transmission coefficient as in area before resonant frequencies (cj < cj^), 
and in the field of resonant frequencies as well. The coefficient of reflecion 
in the field of resonant frequencies has one minimum for thin films (for 
example, for d = 1 nm and v = O.OOlcjp sec^^). 

On Fig. 13 we will show dependence of coefficients of transmission, 
refiection and absorption on quantity of a thickness of a film (i, 1 nm 
^ d ^ 100 nm, hj u = O.OOlcjp sec^^,^ = 75°. Transmittance has one 
minimum in a considered range of thickness, and the absorptance has one 
maximum. 

The formula (24) means, that periodic character of transmittance, refiec- 
tance and absorptance explaints under the presence function tanh(zo/''7o) 
in the second member of this formula. This member is called Debaye mode. 

On Figs. 14 and 15 is shown, that transmittance, refiectance and ab- 
sorptance have extrema in the same points il^, independently of quantity 
of the pitch angle of the electromagnetic wave. These reasons allow to find 
a thickness of a film on those points n„ = oj^/cj^, in which coefficients T, R 
and A have an extremum. 

Let us pass to a deducing of the formula for calculation of the thickness 
of a film in those points il^, in which coefficients of transmission, refiection 
and absorption have extrema. We consider coefficient of refiection. For this 
coefficient on Fig. 16 the first links of the comb represented earlier on Figs. 
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8-10 are considered. In this figure the dot curve corresponding to discrete 
and continuous spectrum, coincides with continuous curve, answering to 
discrete spectrum. 

Points Qn in which the reflection coeflicient has minimum, in accuracy 
coincide with points, in which function cos(Re izo/rjo) possesses the value 
of zero (see Fig. 16). From the equation cos(Re izo/rjo) = we flnd: 

r 77^ r =o+^^' n = 0, 1,2,3, ••• , 

V 77o(^n,£) ^ 2 

or, in exphcit form, 

^ /^ c.,.10-^.(e-in„) A . „ = o,1.2.3,.... (27) 

In the formula (27) the quantity d (thickness of a film) is measured in 
nanometers. The exact formula for calculation of the thickness of a film 
on frequencies Q^, n = 0, 1, 2, • • • is deduced from the formula (27). In the 
points Qn, n = 0, 1, 2, • • • the refiection coefficient has local minima: 

10^7r^;i7(l + 2n) . . ^ . 

d= P J n = 0,l,2,3,---. (28) 

^pRe — -1 

From Fig. 16 it is visible that at n = 3 = 1.025. The formula (28) 
gives the thickness: d = 9.968 nm, i.e. an error in visual determination of 
the thickness of the film gives 0.3%. 

7. Conclusion 

In the present work for the thin films which thickness does not sur- 
pass thickness skin layer, the formulas for calculation of transmittance, 
refiectance and absorptance are received. The analysis of these coefficients 
is carried out. The formula for a finding of a thickness of a film by reso- 
nances in the field of resonant frequencies is deduced. 
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Figure 11: Transmittance, d = 2 nm, u = O.OOScUp. Curves 1,2,3 correspond to values of 
parameters cu = 0.9ujp,uj = ujp,uj = 1.01a;p. 




O.O 0.7 1.4 
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Figure 12: Reflectance, d — 1 nm, u — 0.001a;p. Curves 1,2,3,4 correspond to values of 
parameters cu — 0.95ujp,uj — ujp,u! — l.lujp,uj — 1.2uip. 




I I 

1 10 100 d 

Figure 13: Dependence of transmittance, reflectance and absorptance from fllm thickness 
d, 1 nm ^ d ^ 100 nm, ui — oUp, v — O.OOlaJp. 
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1.0 1.2 1.4 eo/co 



Figure 14: Transmittance, d = 1 nm, v = O.OltUp. Curves 1, 2, 3 correspond to values of 
angle = 75°, 60°, 45°. 
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Figure 15: Absorptance, d — 1 nm, v — 0.01a;p. Curves 1, 2, 3 correspond to values of 
angle 9 = 75°, 60°, 45°. 
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Figure 16: Reflectance, d — nm, v — O.OOlo'p. Periodic curve is set by equation 
y = 0.8 + 0.25cos(Re (i^o/?7o))- 
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